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1. Introduction 

The use of Zernike polynomials for representing the aberrations of an optical system is now 
common and well understood [1]. Indeed, these polynomials have been proven to form a very 
convenient basis in applications as diverse as the description of pupil aberrations of the hu- 
man eye [2] or the characterization in observational astronomy of the statistics of turbulent 
aberrations produced by the atmosphere over the telescope aperture [3]. When using the 
Zernike basis, the aberrated phase is described as a linear combination of the polynomials, 
the value of the coefficients being related to the aperture on which the phase is defined. 
When this aperture undergoes geometrical transformations such as scaling, translation and 
rotation, the Zernike coefficients corresponding to the transformed aperture can be computed 
from proper linear combinations of the original ones by finding the transformation matrix 
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associated to these linear applications. 

In that matter, if many solutions have been proposed for pure pupil scaling, see e.g. [4-8], only 
a few authors have tackled the three linear transformations altogether in a unified method. 
Let us mention Lundstrom & Unsbo [9] who, by generalizing the paper of Campbell [6], 
have developed an ingenious technique that however necessitates a specific reordering of the 
polynomials in the complex plane. Alternatively Bara et al. [10] have proposed a method 
based on the theoretically simple principle of reference change matrix built from a given grid 
of points that sample the wavefront. In practice however this technique appears sometimes 
difficult to implement, especially when the scaling factor is important and/or when the num- 
ber of Zernike coefficients needed to represent the phase is large (the matrix being in these 
cases often ill-conditioned). 

In this paper we propose an alternative technique based on the Fourier Transform properties 
of the Zernike polynomials. We provide analytical and exact formulas that allow the compu- 
tation of the transformation matrix in all cases of scaling/translation/rotation of the pupil. 
In the hypothesis of a full aperture without central obstruction, we develop the results in 
terms of integrals of a product of two or three Bessel functions. The numerical evaluation of 
these integrals can be fastly performed by making use of Mellin Transform properties [11] or 
merely by referring to published tables of integrals [12,13] and it will not be addressed here. 
Section (2) recalls the properties of the Zernike polynomials useful for the formal derivation 
of our calculations. Section (3)provides the general formalism that enablesthe building of the 
transformation matrix, whereas Sects (4, 5, 6), are focusing in further details on the scaling, 
translation and rotation processes, respectively. 

2. Zernike polynomials basic properties 

Any wavefront <j)(r) defined over a circular aperture can be described as a linear combination 
of Zernike polynomials Zj(p): 



where R is the physical radius of the pupil. Zernike polynomials are defined over the circle 
of unit radius n p (p) , that is: 



The Zernike coefficients can be calculated by projecting the phase on the Zernike basis: 




(1) 




(2) 




(3) 
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In polar coordinates, Zernike polynomials are defined for a circular aperture without ob- 
struction as: 



\/2cos(|ra|0) if m > 



Z^(p,9) = Z j (p,9) = Vn + lK(p){ V2sm(\m\e) if m < (4) 

1 if m = 

where n and m are respectively the radial degree and the azimuthal frequency of the j th 
polynomial (note that n and m have necessarily the same parity), j being defined as j = 
n(n + 2)+m + 1, and: 

(n "J^! )/2 f_iW n _,V 

H m (n) = V ' ' [ >' D n ~ 2s (&) 

' ^ s\[(n+\m\)/2-s]\[(n-\m\)/2-s}r {) 

In the full aperture case, Zernike polynomials are orthonormal: 

J U p (p)Z l (p)Z J (p)d 2 p = 5 tj (6) 

and the Fourier Transform Qj(n,a) of U p (p)Zj(p) writes in polar coordinates: 

(-l^n-HV^M^cosflmla) if m > 
Qj(K,a) = (-l) ra v^^ \ (-l)( n "H)/2 z Hy2sin(|m|a) if m < (7) 

■l) n / 2 ifm = 



3. Changing axis reference: building the transformation matrix 

Figure (1) sketches the principle of reference system shifting. The input phase 4>(Rp) is 
defined over the pupil diameter R, in the Oxy reference axis. The output phase (p(R'p') 
defined over the pupil diameter R' is a scaled/translated/rotated version of the input one, 
and is centered on the O'x'y' reference axis. We call t the unitary translation vector between 
both reference systems, and 6 r the rotation angle. For a given point P of the phase that 
belongs to both input and output pupils, we have <f>(P) = 4>(Rp) = 4>(R'p') where: 

Rip' = K 0r [Rp - Rt] ^p' = U 6r [/3(p - 1)] (8) 

with j3 = R/R' and TZe r is the rotation operator. In the following we impose the transformed 
pupil to lie inside the input one (no extrapolation), that is R, R' and t have to verify this 
additional relationship: 

R\t\ +R' < R^ p\t\ + 1 < (9) 
In this framework the input phase is described as a sum of N z Zernike polynomials: 

<f>(R P ) = X> z ;(p) ( 10 ) 




Fig. 1. Sketch of the reference systems used in Sect. (3). In the Oxy reference 
axis (solid lines), the phase is defined over the aperture of radius R and the 
unitary coordinates of a point P in this system is p. In the O'x'y' reference 
axis (dashed lines) the phase is defined over the aperture of radius R' and the 
unitary coordinates of a point P in this system is p' . The reference change 
involves a translation of vector Rt and a rotation of angle 9 r . 



Ragazzoni et al. [14] and Shu et al. [7] have demonstrated that it exists a set of Zernike coef- 
ficients limited by the same highest polynomial number N z that can represent the wavefront 
over the transformed pupil, that is: 

N z 

<P(R>p') = ^Z l (p') (11) 

i=l 

with 

b t = J Il p tf)Z i (p')<f>(R'p')d 2 p' (12) 
Inside the output pupil, we have <ft(R'p') = <p(Rp) and we can modify previous equation as: 

r n z 

h = / n p (p')^(p')0(^P)dV = X>i / ^piff)Zi{ff)Z j {pW (13) 

J j=l J 

that we can rewrite: 

h = 7T % / n p ( P ')^(p')n P (p)^(p)dV (14) 

.7 = 1 J 
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providing that Eq. (9) is verified. Then using Eq. (8) it comes: 

n p (p')^(p')n P (5=|M + 1 ) + t) dV (15) 

Making use of convolution and Fourier Transform properties [15], the latter equation writes: 
ft. = (3 2 n a j / QK K )Qf~ 9r (P K ) ex P [2™Pt.K] d 2 K (16) 

7=1 ^ 

where Qj ° r (K) denotes the Fourier Transform of H p Zj {TZe r [p\). Note that in polar coor- 
dinates IlpZj (JZo r [p\) = HpZj(p,9 — 9 r ), hence its Fourier Transform simply translates as 
Qj 8r (ft, a) = Qj(rz, a — 6 r ). Rewriting Eq. (16) as a product of matrices, we finally have: 

& = M^''Ha (17) 

where a = [ai, ...,a Nz ], b = \b x , ...,6jv.J, and the coefficients (M^' i ' 6 ' r ') of the matrix M^ ,t,e ^ 
verify: 

/•oo /»27r 

M \PtM = n p2 / / «Q*( Kja )Q j .( / g« ja + r ) eX p[2i 7 r / g|t|«cos(a-a t )]dad« (18) 
Jo Jo 

with at = arg(t). Eq. (18) is the general expression of the linear application that enables 
the transformation of the Zernike coefficients in cases of scaling/translation/rotation of the 
pupil over which the wavefront is originally defined. In the following, we will focus on the 
full aperture case without central obstruction such that Eq. (7) can be used for our further 
derivations. 



4. Concentric pupils 



We first investigate the case where the output pupil is a shrunk version of the input one with 
a scaling factor of R'/R =1/(3 (and with \t\ = 0, a t = 0, 9 r = 0). Introducing the expression 
of Zernike Fourier Transform in Eq. (18), it comes: 



M, 







7T 



-(\mj\ + \mj\) 



{\mj\-\mi\) 



y/rii + ly/rij + 1 



x / k 1 J th+1 (2ttk) J nj+1 (2ir/3K)dK 



X 



2tt 



y/2 cos(|mj|a) 
y/2 sin(lmjla) 



a 



\/2cos(| 
y/2 sin(|mj|a) 
1 



da 



(19) 



where the operator (^) denotes all the possible multiplicative combinations between the coef- 
ficients of the vectors, depending on the respective values of m; and rrij (i.e. 9 multiplications 
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Fig. 2. Center: initial phase (top) and associated Zernike coefficients (bottom), 
with j = 1 to 28, that is the phase is defined with n = 6 radial degree modes. 
Left: reconstructed phase over the scaled aperture (/3 = 1.3) computed from 
the transformed Zernike coefficients. These coefficients are calculated from the 
original ones using both the transformation matrix of Eq. (20) (circles) and 
Shu et al. method [7] (solid line) . Right: same as previously for a scaled and 
translated aperture ((3 = 1.5, \t\ = 0.3, a t = tt/3). The transformed coefficients 
derived from Eq. (23) (circles) are compared with that of Lundstrom & Unsbo 
technique [9] (solid line). 



possible in this case). From the trigonometric properties of cosine functions, it comes straight- 
forward that the integral over a in Eq. (19) is non zero only if Wj = rrij. Hence previous 
equation can be replaced by: 



2f3y/rii + ly/rij + 1 / °° k 1 J nt+ i(27iK)J n . + i(2'rr(3K,)dK if m 



m 



if Mj 7^ m 



j 
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Figure (2, left and center) shows an example of pupil scaling with /3 = 1.3. We can verify on 
the upper plots that the computed phase over the scaled pupil is indeed the zoomed version 
of the input one, with the appropriate zooming factor /3. Bottom graphics are displaying the 
input and derived output Zernike coefficients. In order to further check the validity of our 
method, we have also plotted the estimated output coefficients calculated with that of Shu 
et al. [7]. Both independent techniques are giving the same results. 

5. Scaling and translating 

Since we have assumed that the transformed pupil should lie within the original one, any 
shifting must necessarily be combined with appropriate scaling, such that Eq. (9) is verified. 
As a consequence, pure translation process cannot be tackled by this study unless a so-called 
meto-pupil encompassing both pre- and post-shifted apertures is considered beforehand [14] . 
Providing that this hypothesis is verified, Eq. (18) rewrites (with 9 r = 0): 

r \P,t] P f ■Mn i +n J )-\rr H \- 



M. 



'j 



7T 



2n 



\/2 cos(|raj|a:) 
v2 sin(|mi|a) 
1 



i(KI-KI)V^nV^+ T / ^ 1 ^ +1 (27TK)J ni+1 (27r/3« ; ) 

Jo 

exp [2z7r/3|t|ft cos(o! — a t )} dadn (21) 



a/2 cos(|mj|a) 
a/2 sin (| rrij \ a) 



The derivation of the integral over a requires to make use of the integral definition of Bessel 
functions [16], that is: 

1 /"7T 

(22) 



J n ( K ) = — / ex P [~ i(na — Ksina)] da (with n integer) 

The formal development of Eq. (21) necessitates several steps that are not difficult however 
somewhat tedious. They are presented in App. (A). The coefficients of the matrix M^'*! 
eventually write: 

= V2° 



M ■ — a/2 3 ' (—\\ — — 



X 



A 



-A + 



k 1 J ni+ i(2nK)J ni+1 (27if3K)J\\ mz \_\ nh \\(2Tif3\t\K)dK 
K _1 J ni+ i(27rK)J nj+1 (27r/3K)J| fni | + | mj |(27r/3|t|fi:)]o?fi: 



(23) 



with 



A- 



and 



A 



llmjl-lmjii J cos([|mj| — |mj|]o;t) if sgn(mj) = sgn(rrij) 

sgn(mj) sin([|mj| — |mj|]a t ) if sgn(mj) ^ sga(rrij) 



I sgn(mj) cos([|mj| + |mj|]at) if sgn(mj) = sga(m,j] 
sin([|mj| + |mj|]a t ) if sgn(mj) ^ sgn(m,j] 
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Input phase Rotated pupil 




Normalized radius Normalized radius # Zernike polynomials 



Fig. 3. Left: original wavefront. Center: reconstructed wavefront for a pupil 
rotation of 9 r = ir/3. Right: transformed Zernike coefficients derived from the 
original ones using Eq. (28) (circles) and Bara et al. method [10] (solid lines). 



If Eq. (23) appears to return complex numbers, the conjugated parity of rij and m, (respec- 
tively rij and rrij) insures to produce a real coefficient for any pair of Zernike polynomials 
considered. Figure (2, center and right) illustrates the case of scaling + translation with 
(3 = 1.5, \t\ = 0.3, and 9 t = 7r/3 and results of the present method are proven to match with 
that of Lundstrom & Unsbo [9] . 



6. Rotation 

We finally investigate the case of pure rotation as it can be treated independently from the 
scaling and translation processes. A full scaling+translation+rotation transformation can 
indeed be decomposed in two consecutive processes, that is scaling+translation then pure 
rotation (or similarly rotation first, then scaling+translation). In other words the transfor- 
mation matrix of the full transformation M^'*' 6 ^ can be computed as: 



M [/3,tA] _ mWA.M^ = M [6r] .M [/3 ' t] 
Applying Eq. (18) with /3 = 1 and \t\ = 0, it comes: 



M. 



-(\m i \ + \m A) 



'J 



1 3(»j+»,-) 

- (-1)— 3— -1 s z^H^V^+V^Tl 

7T V 



X / K 1 J ni+ i{2'KK)J n +1 (27TK)dK 



2tt 



a/2 cos (| mj| a) 
y/2 sin (| rrii \ a) 



a/2 cos ( | rrij \ [a + 9 r \) 
y/2s,m{\mj\[a + d r \) 



da 



(24) 



(25) 
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In this case, trigonometric properties of cosine functions imply that the integral over a is 
non-zero only when \rrii\ = \rrij\, and Eq. (25) reduces to: 



2y/ui + lyjrij + 1 / k 1 J ni+1 (2ttk) J n j+1 (2'KK)dn 
Jo 



cos(rrij6 r ) if 



sin(rrij9 r ) if = — rrij (and 7^ 0) 
if \rrii\ 7^ \m,j\ 



(26) 



Furthermore it can be shown that: 



POO 

2y/m + l^Tlj + 1 / K~ l J ni+ i(2TXK) J nj+1 {2-KK)dK = 5, 

Jo 

so that M\j takes the simpler form: 

if m 



n,i,rij 



(27) 



*3 



cos(rrij9 r ) 
sm(mj6 r ) if 



-mj (and m 8 7^ 0) 



(28) 



if \rrii\ ^ \m,j\ 

As a consequence, the rotation matrix does not necessitate the evaluation of integral of Bessel 
functions but merely requires to compute cosine functions for multiples of the rotation angle 
accordingly to the accounted azimuthal frequencies (rrij). Figure (28) shows an example of 
phase rotation with 6 r = 7r/3. The transformed Zernike coefficients calculated from Eq. (28) 
are successfully compared with that of Bara et al. [10]. 

7. Summary 

This paper provides analytical formulas that enable the computation of Zernike coefficients 
for scaling, translation and rotation of the aperture over which the wavefront is originally 
defined. The expressions of the transformation matrix that allows to switch from the original 
Zernike coefficients to the transformed ones are exact. Each coefficient of the matrix can be 
calculated directly and independently without involving any specific ordering of the poly- 
nomials. Our method has been compared to several independent techniques avalaible in the 
literature. The results are consistently matching, hence validating the approach presented in 
this paper. 

Appendix A: Computing the scaling+translation matrix 

We focus here on the integral over a of Eq. (21): 



F(k) 



2- 



y/2 cos(|mj|a) 
\/2 sm(\m,i\a) 
1 



a/2 cos (I rrij \ a) 
v2 sin (| rrij | a) 



exp [2i7r/3|t|« cos(a — a t )] da (A-1) 
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We recall the integral definitions of Bessel functions that will be used to unfold Eq. (A-l): 

f 2n / \ i ■ i w a \ 2vr(-l)^ cos(ma t ) J\ m \(y) if m even 
/ cos{ma) exp(ty cos{a - a t ))da = < ^h-i ' (A-2) 

Jo 2z7r(— 1) 2 cos ma t J\m\{y) it m odd 



f 2 * ■ t \ c i \\a / 2tt(-1) I 2 I sin(mai)J H (?/) if m even 

Jo 2z7T — 1 2 sin maiUui « it m odd 



A-27T 
< 
o 2z7r(— 1) J ^ — sin(maj) J| m |(y) if m odd' 

We now have to investigate the different combinations regarding mj and m^. 

1. mj > and m^- > 0: Eq. (A-l) rewrites: 



F(k) = I cos([\mi\ — \mj\]a) exp[2iir/3\t\Kcos(a — a t )]da 
Jo 

+ / cos([\mi\ + \mj\]a) exp[2i7r(3\t\K,cos(a — a t )] da 
Jo 



that, according to Eq. (A-2) gives: 
(a) if rrii and mj have same parity: 



| — \m.j | 

F(k) = 2tt(-1) a cos([|77i<| - |m 3 -|]a t )J|| m( |_| mi ||(27r^|*|«) 

| \m i \ + \rrij \ \ 

+ 2vt(— 1) 2 cos([|mi| + |mj|]o£t) J\\ mi\ + \mj ll(27r^|t| K ) (A-4) 
(b) if mj and m^ have different parity: 



F(k) = 2i7r(— 1) 5 cos([|mi| - |mj|]at)J|| mi |_| TOi ||(27r^|t|K) 

| |mj | + jmj | | — 1 

+ 2z'7r(-l) a cos([|mj| + |mj|]at)J|| TOi | + | mj ||(27r^|t|/c) (A-5) 
2. mj < and mj < 0: 

/■27T 

F(k) = / cos([|mj| — |raj|]a) exp [2i7rf3\t\KCOs(a — at)] da 
Jo 



2tt 



cos([|mj| + |mj|]a) exp [2i7r/3|i|ftcos(a — a t )] da 



o 



that rewrites: 



(a) if mj and mj have same parity: 



\m i \-\m i 



F(k) = 2tt(-1) 2 cos([|mj| - |m 3 -|]a t ) J|| mi |_| mi ||(27r/3|*| 



N™j|+|™jN 

2tt(-1) 2 cos([|mj| + |m i |]a t )J|| mi | + | mi ||(27r/?|t|K) (A-6) 
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(b) if rrii and rrij have different parity: 

F(k) = 2z'7r(-l) 2 cos([|mj| - \mj\]a t )J\\ mi \-\ mj \\(2ir/3\t\K) 

- 2z'7r(-l) 2 cos([|mt| + |mj|]at)J|| m< | + | mj .||(27r^|t|Ac) (A- 7) 

rrii < and > 0: 

F(k) — — I sm([\rrii\ — \rrij\]a) exp[2iir(3\t\K,cos(a — a t )} da 
Jo 

+ / sin([|mj| + |mj|]a) exp [2m/3|t|/t cos(a — a t )] da 
Jo 

that, according to Eq. (A-3) gives: 

(a) if rrii and rrij have same parity: 

l m iM m jl 

F(k) = -2tt(-1) a sin([|mi| - |mj|]at)J|| m4 |_| mi ||(27r^|t|/c) 

+ 2tt(-1) 3 sin([|mi| + |m i |]Q; t )J|| mi | + | m .||(27r / S|t|K) (A-8) 

(b) if rrii and rrij have different parity: 

| |mj | — |?7ij | | — 1 

= -2m(-f) a sin([|mi| - |m i |]o; t )J|| mi |_| mi ||(27r/?|t|«;) 

+ 2nr(-l) 2 Bin([|m i | + |m i |]a!t)J|| m4 | + | mi ||(27r^|t|K) (A-9) 

rrii > and rrij < 0: 

/.2tt 

-^X^) = / sm({\rrii\ — \rrij\]a)exp[2i7r(3\t\Kcos(a — a t )]da 
Jo 

+ / sin([\rrii\ + \rrij\]a) exp[2iir/3\t\K,cos(a — a t )] da 
Jo 

(a) if rrii and rrij have same parity: 

| | — \m,j \ \ 

F(k) = 2tt(-1) 3 sin([|mi| - |m y |]a t )J|| m .|_| mi ||(27r/3|t|«;) 

ll m il + l m jll 

+ 2vr(-l) 2 sm([|mi| + |m i |]a! t )J|| m .| + | mj ||(27r / 9|*|«) (A-10) 

(b) if rrii and rrij have different parity: 

F(k) = 2i7r(— 1) 2 sin^flmil - |m_j-|]a; t )J|| mi |_| mj .||(27r^|t|K) 

+ 2ztt(-1) 2 sin([|mi| + |mj|]a t )J|| m .| + | mi ||(27r/3|t|/«) (A-ll) 
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5. rrii > and rrij = 0: 



1>Z1T 

F(k) — \[2 / cos(|mj|a) exp [2wr/3|t|« cos(a — at)] da 
Jo 

(a) if rrii is even: 

F{k) = 2nV2(-l)— cos(|m i |at)J K |(27r/3|t|«) (A-12) 

(b) if rrii is odd: 

F(k) = 2W2(-1)^ cos(|m<|at) J| fn .|(27r^|<|«) (A-13) 
6. rrii < and = 0: 



F(k) — v2 / sin(|mj|a) exp [2z7r/3|i|re cos(a — a t )] <ia 

(a) if ?Tij is even: 

F(k) = 2ttv / 2(-1)— sin(K|a t ) J\ ml {2Trf3\t\K) (A-14) 

(b) if rrii is odd: 

F(re) = 2mv / 2(-l)^ sin(|m i |at)J K |(27r/3|*|K) (A-15) 

7. mi = and rrij > 0: 

F(k) = V2 / cos(|mj|a) exp [2z7r/3|i|re cos(a — a t )\ da 
Jo 

(a) if rrij is even: 

F(re) = 2 7 rv / 2(-l)— cos( I rrij | a t ) J K , (2nf3\t | re) (A-16) 

(b) if 77ij is odd: 

F(re) = 2my/2{-VT^ 1 cos(|m,-|at) J| m .|(27r/3|*|re) (A-17) 

8. rrii — and rrij < 0: 

F(re) = v2 / sin(|mj|a) exp [2i7r/3|t|re cos(a — at)] da 



(a) if rrij is even: 



|m ■ I 



F(re) = 2ttv / 2(-1)^ 1 " sin(|m y |at) J K |(27r/3|t|re) (A-18) 
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(b) if rrij is odd: 

|mJ-l 

F(k) = 2ctv / 2(-1)^ _ Bindmjlot) J| mj .|(27r/3|*|«) (A-19) 
9. vrti = and rrij = 0: 

/■2tt 

F(k)= exp[2i7r/3\t\Kcos(a- a t )) da = 2nJ (2ix (5\t\K) (A-20) 

Putting all together the different combinations leads to the general expression of Eq. (23). 
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